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The dynamics of the domains is studied in a two-dimensional model of the microphase separation 
of diblock copolymers in the vicinity of the transition. A criterion for the validity of the mean field 
theory is derived. It is shown that at certain temperatures the ordered hexagonal pattern becomes 
unstable with respect to the two types of instabilities: the radially-nonsymmetric distortions of the 
domains and the repumping of the order parameter between the neighbors. Both these instabilities 
may lead to the transformation of the regular hexagonal pattern into a disordered pattern. 
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Formation of complex patterns consisting of do- 
mains with sharp walls is a beautiful example of self- 
organization in the systems both near and far from 
thermal equilibrium Recently, it became clear 

that long-range competing interactions are responsible 
for the formation of the domain patterns in the systems 
as diverse as ferroelectrics and ferrofluids, garnet ferro- 
magnets, Langmuir monolayers, type I superconductors 
in the intermediate state, diblock copolymers, reaction- 
diffusion systems with long-range inhibition (see Q and 
references therein) . In such systems the formation of the 
uniform state favored by the local properties of the sys- 
tem is precluded by the long-range interaction which does 
not favor that uniform state. Thus, the system becomes 
separated into the domains of the alternating values of 
the "order parameter" . As a result of this separation a lot 
of different equilibrium configurations, including highly 
symmetric ones, are possible. 

In this Letter we will study the time-dependent model 
of a system with competing interactions in two dimen- 
sions. We will investigate the stability of the station- 
ary states and show that they undergo instabilities which 
may change both the characteristic length scale and the 
morphology of the domain patterns. We will also show 
that the destabilization of the highly symmetric patterns 
typically leads to the formation of the disordered pat- 
terns. 

The microphase separating diblock copolymer melts 
are a typical example of the systems with the long-range 
competing interactions. There the macroscopic phase 
separation of the mutually incompatible monomers is not 
allowed since the monomers are connected through the 
polymer chains. Ohta and Kawasaki obtained the free 
energy for this system in the case of the long polymer 
molecules 
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+ ^ / dV(0(x) - <t>)G{x ~ x'mx') - 0) , (1) 



where 4> is the order parameter which is proportional to 
the difference of the monomer concentrations, a and h 
are the coefficients of the Landau-Ginzburg expansion, 
the last term in Eq. (^ is the long-range interaction 
characterized by the function G which reflects the con- 
nectivity of the chains, a ~ N~'^ is the strength of this 
interaction and N is the number of monomers in a chain, 
the constant cj) is determined by the block ratio. The 
function G satisfies 
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The model given by Eqs. (|l]) and (^) in fact has a wider 
applicability and can be used to describe the station- 
ary states in ceramic compounds with the long-range 
Coulombic interactions ferroelectric semiconductors 
Q, high-temperature superconductors and degenerate 
magnetic semiconductors ||T^, reaction-diffusion sys- 
tems with the long-range inhibitor pjlll ] , and reaction- 
controlled spinodal decomposition 

The formation of the domains in the system under con- 
sideration is due to the competition between the nonlocal 
interaction and the surface tension which is determined 
by the local terms in Eq. (|l|). For the equilibrium pattern 
the contributions of these two effects have to be compa- 
rable. According to Eq. (|l|), for the domain of size R we 
have 
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where a is the coefficient of the surface tension. Near the 
critical point a = (To\t\''^'^-'^\ where t = {T-Tc)/Tc, is 
the critical temperature, and v is the critical exponent of 
the Ising model jl^ . Because of the long-range character 
of the nonlocal interaction, in estimating its contribution 
one can ignore the fluctuations of (f> and put (j) — (t>o\t\^ 
in Eq. (m, where (3 is the respective critical exponent of 
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the Ising model. Then, according to Eq. the char- 
acteristic size of the equihbriiim domain near the critical 
point will be 
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Substituting this expression into Eq. , we obtain that 
the free energy of a single domain AF ^ 1 , if 
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If this condition is satisfied, we can use the mean-field 
approximation to describe the domain patterns in the 
vicinity of the phase transition. 

In order to study the domain patterns in the mean-field 
limit let us introduce the new variables 
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where ^ — Co 1^1 " is the correlation length, and write 
the free energy functional (up to a constant factor and 
dropping the primes) 

+ ^ j d'x'Wx) ~ ^)G{x - x'mx') - 0)^ . (7) 

One should think of Eq. (0) as the mean-field represen- 
tation of the domain interactions. When the condition 
of Eq. (H) is satisfied, we have e ^ 1, what implies the 
strong segregation limit 0. Notice that in these units 

The evolution of the oder parameter (j) near the phase 
transition in the mean-field limit can be described by 
the time-dependent Landau-Ginzburg equation obtained 
from Eq. (|^) . This equation can be written as a pair 
of reaction-diffusion equations of the activator-inhibitor 
type in the limit of the fast inhibitor: 
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where (j) plays the role of the activator and V' plays the 
role of the inhibitor and the time scale has been ab- 
sorbed into the definition of t. In the limit e ^ these 
equations are equivalent to the interfacial dynamics prob- 
lem pT] , |T5| . We would like to emphasize that in this 
limit Eqs. (||) and (j^) are applicable to the systems with 
the conserved order parameter. Indeed, the equations of 
the interfacial dynamics mentioned above can be alter- 
natively derived from the free energy given by Eq. if 
one assumes the overdamped dynamics of the interface 
[ pT[ . The latter is justified for the microphase separation 
in the diblock copolymer melts, where the hydrodynamic 
effects are important. 



The domain patterns that form in this kind of systems 
have been recently studied in a great detail in connection 
with the patterns forming in highly nonequilibrium sys- 
tems [p^-p^. Muratov and Osipov showed that in the 
case of Eqs. (^) and (||) only the patterns whose charac- 
teristic length scale (in these units) is of order e~^/'^ can 
be linearly stable p5|-p7[|. Thus, if the condition in Eq. 
(^ is satisfied, on this length scale the mean-field ap- 
proximation remains valid for the time-dependent theory 
as well. 

When 1^1 < 1 and e < 1, Eqs. (|) and (|) admit 
the solutions in the form of the stationary multidomain 
patterns, including the hexagonal patterns of circular do- 
mains of radius TZs and the period Cp in two dimen- 
sions ^,[7|,p^-[l8| . It is clear that the period of these 
patterns (or, actually, the characteristic interdomain dis- 
tance) must lie within I <^ Cp ^ L, where I ~ 1 and 
L ^ e^^ are the characteristic length scales of the varia- 
tion of (p and ip, respectively [Q. For \(j>\ not very close to 
1 the radius of the domains is comparable to Cp. Since 
the characteristic size of the domains is of order e^^^^, 
the period of the pattern must also be Cp ^ e"^/'^ <C L, 
so the value of (f> is close to -1-1 inside the domains and 
close to —1 outside. 

Let us now study the stability of the hexagonal pat- 
tern. It is clear that the results obtained for this pattern 
should be qualitatively the same for an arbitrary mul- 
tidomain pattern. In this case, according to Eq. (|^), 
the relationship between the period and the radius of the 
domains is 
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so the ratio TZs/ Cp can be conveniently used as a param- 
eter instead of 0. According to the general asymptotic 
theory of instabilities of domain patterns in reaction- 
diffusion systems the dangerous fluctuations of (j) 
are localized in the domain walls and represent the small 
distortions of the domains. The damping decrement 7 
of such fluctuations is determined by the eigenvalues of 
a certain operator. Since the hexagonal pattern we are 
interested in is periodic, we can partially diagonalize this 
operator by considering the fluctuations modulated by 
the wave vector k which lies in the Brillouin zone and 
obtain for the system under consideration |l6| 
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where m and m! are integers corresponding to the modes 
of the azimuthal distortions of the domains, Aq is a con- 
stant independent of m and fc, 5mm' is the Kronecker 
delta, and Rmm'{k) is a certain matrix with the indices m 
and m'. It is possible to show (the details of the deriva- 
tions will be given elsewhere) that the constant Aq is 
given by 
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Since the characteristic length scale of the pattern is 
much smaller than L, the Laplacian dominates in the 
equation for the Green's function involved in the calcula- 
tion of Rmm'{k) (see Ref. so this Green's function 
may be assumed to satisfy Eq. (H). With this fact in 
mind it is convenient to write the expression for Rmm' (k) 
in the Fourier space 
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where k' runs over the reciprocal lattice, Jm are the 
Bessel functions, "dk+k' is the angle between the vector 
k + k' and the x-axis, and v = \/3£p/2 is the volume of 
the elementary cell. Notice that for TZs <C Cp the off- 
diagonal elements of Rmm' {k) become negligibly small, 
so for |m| > 2 we recover the instabilities of the local- 
ized domains — autosolitons |l^]. As TZs gets bigger, the 
mixing between the vectors corresponding to the differ- 
ent values of m occurs. This mixing, however, is not very 
strong, so one can still label the eigenvectors of Eq. ( |ll| ) 
with m. 

As was shown qualitatively by Kerner and Osipov, for 
the most dangerous fluctuations the wave vector k will 
lie close to the edge of the Brillouin zone Q. There are 
two basic types of the fluctuations we need to consider: 
the fluctuations with m = which lead to repumping 
of the order parameter between the neighboring domains 
(Fig. |l|[a]) and the fluctuations with m = 2 which lead to 
the asymmetric distortions of the domains (Fig. |l|[b]) 
The analysis of Eq. ([l^) shows that the most dangerous 
fluctuations with m = have k = i(&i — 62), where bi and 
62 are the reciprocal lattice vectors which make a 120° 
angle, while the most dangerous fluctuations with m = 2 
have k — ^(^1 + ^2)- The instability 7 < occurs with 
respect to repumping when Cp < Cpo or with respect to 
the asymmetric distortion when Cp > Cp2, where >Cpo,2 
depend on e and TZsjCp. The resulting stability diagram 
is presented in Fig. Figure || also shows the period of 
the equilibrium hexagonal pattern obtained by Ohta and 
Kawasaki 0. One can see that the equilibrium hexago- 
nal pattern is stable for all values of TZs (except, possibly, 
for TZs/Cp close to 0.5 where the assumption about the 
circular shape of the domains ceases to be valid). In 
two dimensions the equilibrium hexagonal pattern cor- 
responds to the global minimum of the free energy if 
Hs/Cp < 0.31, whereas for 0.3_1 < TZs/Cp < 0.37 (the 
second condition means that < 0, i.e., positive do- 
mains in the negative background) the global minimum 
corresponds to the lamellar pattern [Q. Figure ^ how- 
ever, does not show the transition from the hexagonal to 



the lamellar pattern, so in fact the equilibrium hexagonal 
pattern is always metastable. 

Recall that the values of e and (f> strongly depend on 
temperature near the critical point. Suppose that the 
equilibrium hexagonal pattern formed as a result of the 
slow quench of the system below Tc- If now the temper- 
ature is abruptly raised, the domains may become un- 
stable with respect to the asymmetric distortions. To 
study the kinetics of this process we solved Eqs. (^) and 
(^) numerically with the initial condition in the form of 
the hexagonal pattern plus small noise. This process is 
shown in Fig. ^(a). At the end of the simulation the 
system reaches an asymptotic state. One can see that 
the hexagonal pattern transform into a highly disordered 
metastable pattern which consists of the domains of com- 
plex shapes. Notice that this effect was observed experi- 
mentally in different systems with competing interactions 

If the system temperature is quickly lowered, the pat- 
tern may destabilize with respect to repumping. The ki- 
netics of this process is shown in Fig. §(b). As the desta- 
bilization progresses, some of the domains start to "eat" 
their neighbors, what results in either increase of the in- 
terdomain distance or fusion of some of the neighbors. A 
lot of the domains shrink and eventually disappear. All 
these processes create a lot of disorder. Eventually the 
distance between the domains becomes large enough, so 
that the repumping is no longer realized, and the result- 
ing pattern orders somewhat at late stages. However, 
the asymptotic pattern in the end of the simulation is 
still highly disordered. Notice that the repumping effect 
is observed in experiments and numerical simulations of 
systems with competing interactions |^,^ . 

In summary, we showed that the stationary multido- 
main patterns in the systems with long-range compet- 
ing interactions may undergo instabilities affecting both 
the period and the shape of the domains which typically 
result in the formation of highly disordered metastable 
patterns. 
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FIG. 1. Two major types of instabilities of the hexagonal 
pattern. 

FIG. 2. Stability diagram for the hexagonal pattern. The 
dashed line corresponds to the equilibrium pattern. 



FIG. 3. (a) Destabilization of the hexagonal pattern with 
respect to the asymmetric distortions. Distributions of 
at different times for e = 0.05 and = 0. The system is 
200 X 230. (b) Destabilization of the hexagonal pattern with 
respect to repumping. Distributions of 4> at different times 
for £ = 0.025 and = -0.2. The system is 200 x 230. 
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